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ABSTRACT 

We present results from a numerical forward model to evaluate one-dimensional reduced power 
spectral densities (PSD) from arbitrary energy distributions in k-space. In this model, we can 
separately calculate the diagonal elements of the spectral tensor for incompressible axisymmetric 
turbulence with vanishing helicity. Given a critically balanced turbulent cascade with k\\ ~ 
and a < 1, we explore the implications on the reduced PSD as a function of frequency. The 
spectra are obtained under the assumption of Taylor’s hypothesis. We further investigate the 
functional dependence of the spectral index k on the field-to-flow angle 9 between plasma flow 
and background magnetic held from MHD to electron kinetic scales. We show that critically 
balanced turbulence asymptotically develops toward ^-independent spectra with a slope corre¬ 
sponding to the perpendicular cascade. This occurs at a transition frequency / 2 d(T, a, 9), which 
is analytically estimated and depends on outer scale L, critical balance exponent a and field-to- 
flow angle 9. We discuss anisotropic damping terms acting on the k-space distribution of energy 
and their effec t s on t he P SD. Further we sh ow that the spectral anisotropies k(0) as found by 
Horburv et al. (12008 1 and Chen et ahl ( 20101 1 in the solar wind are in accordance with a damped 
critic ally balance d casc ade of kinetic Alfven waves. We also model power spectra obtained by 
von Papen et ah (|2014h in Saturn’s plasma sheet and find that the change of spectral indices 
inside 9 R s can be explained by damping on electron scales. 


Subject headings: turbulence — solar wind — plasmas — planets and satellites: magnetic fields — 
methods: numerical 


1. Introduction 


Plasma turbulence has been analyzed exten¬ 
sively using in-situ measurements of magnetic 
field and velocity fluctu a tions i n the sola r wind 


(Matthaeus & Goldstein 1982; Tu et al. 

1984 

Burlaea & Klein 1986; Bruno & Carbone 

200511. 


The statistical properties of these fluctuations 
are usually interpreted by means of power spec¬ 
tral density P(f) as a function of frequency in 
spacecraft frame. Power spectral densities (PSD) 
derived from these in-situ measurements of sin¬ 
gle spacecraft are generally obtained in a one¬ 
dimensional reduced form. This means that fluc¬ 
tuations associated with various wave vectors are 
observed at the same frequency in spacecraft 
frame and thus contribute to the spectral den¬ 


sity P(f) at a certain frequency /. However, 
turbulence models are commonly formulated in 
three-dimensional wave vector space, where each 
wave vector k is assigned a certain energy. Under 
the assumption of Taylor’s hypothesis and sta¬ 
tistical homogeneity, the observable power spec¬ 
trum P(f) can be obtained from the energy dis¬ 
tribution in k-space by a two-dimensional integra- 
tion. Although this f act is very well known (e.g . 


Fredricks fe Coroniti 1976: Forman et al. 201 it 


Wicks et, al.ll2012 ; Turner et, al. 2012lh there is cur¬ 


rently no application, where several plasma scales 
in k-space have been considered simultaneously. 
We provide here a numerical tool, which we use to 
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study the effects of various k-space distributions Q 

Note, that for the calculation of reduced power 
spectral densities in frequency space, the appli¬ 
cability of Taylor’s hypothesis and homogeneity 
needs to be justified for each physical system un¬ 
der investigation, e.g., the solar wind or magneto- 
spheric plasmas. We will return to this issue when 
we compare observations in the solar wind and 
Saturn’s magnetosphere with our model results. 
The transformation from k-space to the observed 
frequency space taken in a rest frame moving with 
respect to the turbulent media is unique. In con¬ 
trast, the observations taken in the moving frame 
do not allow to uniquely determine the spatio- 
temporal turbul ent str u cture in the rest f rame of 
the media (e.g. Fredricks fe Coronitil 1976h . Thus 
interpretation of observed single spacecraft fre¬ 
quency spectra can only test whether certain k- 
space distributions of turbulence are consistent 
with observations or not, but they cannot uniquely 
identify one certain k-space structure. 

The statistical properties of turbulent fluctu¬ 
ations have led to the understanding that tur¬ 
bulent fluctuations in plasmas are anisotropic 
with respect to a back ground magnetic field B 0 . 
Mat thacus et al. ( 199C)h observed anisotropy with 
a composite slab and 2-D model consisting of wave 
vectors parallel and perpen dicu lar _to the mean 
magnetic field, respectively. Bieber et al. ( 1996i l 


ical balance is an inhere n t pro perty of plasma 


turbulence. iTessein et al. ( 20091) ■ e.g., found no 


successfully explained the relative ratios of turbu¬ 
lence power in different magnetic field components 
of Helios m easurements with a 85 % 2-D and 15% 
slab model. ISaur fe Bieber ( 19991 ) found that low 
frequency turbulent fluctuations in the solar wind 
could be explained with a combination of 2-D and 
a radial slab model. To expl ain the ani so trop ic 


cascading process, Goldreich fe Sridhar ( 19951 ) 


proposed a critical balance between Alfven wave 
period and nonlinear time, which leads to spectral 
anisotropy of the measured power spectral den¬ 
sities. Although there have been recent results 
that seem to confirm the critical balance theory 
(IHorburv et al. 2008 : Podesta 2009t Wicks et al 


2010; C hen_et al.1 2010t TenBarge fe Howes 2012 : 
He__et al- 120131) . it is an ongoing debate if crit- 


1 The numerical code is made available in 
the supplementary online material and under 
www.geomet.uni-koeln.de/en/research/turbulence 
to be used for other k-space distributions by the interested 
reader. 


dependence of the spectral index on the field-to- 
flow angle 0 , i.e., the angle between Bo and plasma 
flow v (relative to the observing spacecraft), using 
a global backgr ound magn e tic fie ld in their analy¬ 
sis an d iGrappin fe Mullen (l2010t) recently found a 
^-independent range, which they claim contradicts 
the critical balance assumption. 

In this paper, we analyze in detail the spec¬ 
tral form of one-dimensional reduced PSD based 
on critical balance in the MHD and kinetic range 
of scales. The PSD are modeled for the first time 
for oblique angles 0° < 8 < 90° from MHD to 
electron scales based on a three-dimensional en¬ 
ergy distribution. In section [2j we discuss the 
diagonal elements of the spectral tensor, which 
depend on two scalar functions for toroidal and 
poloidal fluctuations and we show how the ten¬ 
sor reduces to a one-dimensional spectrum if Tay¬ 
lor’s hypothesis holds. In section [3l we introduce 
our model for the energy distribution in k-space 
based on a critically balanced cascade on MHD, 
ion kinetic, and electron kinetic scales. Further, 
we discuss potential anisotropic damping terms. 
In a comparative study i n section |U we analyz e if 
the results published by Horburv et ahl ( 20081) in 
the MHD regime and Chen et al. ( 2010l) in the ki¬ 
netic regime of the solar wind can be explained by 
a critic ally balanced or slab /2-D model. In addi¬ 
tion to iForman et al. ( 20 111), who made a simila r 
analysis for the results of IHorburv et al.l ( 2008 ). 
we include kinetic scalings and damping effects. 
In section El we investigate how spectra observed 
by von Paoen et al.l ( 20141 ) in Saturn’s magneto¬ 
sphere are consistent with a critically balanced 
cascade and damping processes. 

We like to point out that our analysis is also 
relevant to modeling cosmic ray transport in tur¬ 
bulent media (e.g., Jokipii, 1966; Engelbrecht 
& Burger, 2014, 2015). Cosmic rays are trav¬ 
eling with very large velocities through turbu¬ 
lent plasma fields while being scattered by time- 
variable fields, which are seen in a Lagrangian 
frame by the moving particles commonly described 
in an appropriate guiding center. Thus, cosmic 
rays are subject to the same sampling effects of 
anisotropic wave vector distribution as the time- 
dependent spacecraft measurements commonly 
also taken at large velocities with respect to the 
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turbulent plasma. The effect of field-to-flow an¬ 
gle on cosm ic ray properties h as for example been 
studied by Bieber et al. ( 1996h . 


2. Diagonal Elements of the Spectral Ten¬ 
sor 


As we are primarily interested in the power 
spectra of the fluctuations, we will only dis¬ 
cuss the diagonal elements of the spectral ten¬ 
sor. Under the assumption of incompressible 
axisymmetric and mirror symmetric turbulence 
with regards to the mean magnetic field, a gen¬ 
eral form of the correlation tensor can be derived 


(Matthaeus fe Smlthlll981; Montgomery fe Turner 


1981 ). Mirror symmetry leads to vanishing helic- 
ity (Oughtpn et al. 1199711 and is therefore equiva¬ 
lent to balanced turbulence, i.e., the same amount 
of energy in waves parallel and anti-parallel to the 
background mag netic fi eld (for the imbalanced 
case, see iLithwick fe Goldreichl l2003i l. The di¬ 
agonal elements of the spectral tensor then de¬ 
pend on only two independent scalar functions T 
and <f>, which describe the toroidal (shear-Alfven 
mode on MHD scales) and poloidal fluctuations 
(pseudo Alfven mode), respectively. The latter 
are the incompressible limit of the slow mode 
( Goldreich fe S rid 1 i a.ij 1 1 99 5l ICho et aJ1l2002 1 . 


In a magnetic field-aligned coordinate system 


with Bn 


and e x and e y arbitrarily ori¬ 


ented perpendicular unit vectors perpendicular 
to e . the diagonal elements can be written as 


(jOughton et al.lll997l : IWicks et al.l l2012f 


S xx ( k) 

Syy(k) 

S zz ( k) 


kl 
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k 2 k 2 
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( 2 ) 

(3) 


where k \ = k 2 + k 2 and fc|| = k z . The trace of the 
tensor is tr(5”(k)) = + $. Because the general 

form of U/ and *I> is not known, we need to intro¬ 
duce reasonable estimates for t he scalar fu nctions 
T and *I>. According to iLithwick fe Goldreichl 


(|200ll ), slow mode waves are passively cascaded by 
Alfven waves, which means that the same scaling 
can be expected for the poloidal fluctuations $. 
Due to the fact that k ± 3> fen, which is commonly 


observed in plasma turb ulence (IMatthaeus et al 
Il99f)l : iBieber et al. 1996), we may estimate the 
proportionality between the two functions from 
the measured power anisotropy 


P ± _ S ± _ k 2 T _ T 


(4) 


where P_l = P X x + Pyy is the perpendicular 
and P|| = P zz the parallel component of the 
one-dimensional reduced spectrum. Observations 
in t he solar wind i ndica t e that usually P \ 

P|l ( Hamilton et al. 2008 : TenBarge et al. 2012 ). 
which means the error stemming from the assump¬ 
tion of passively cascading poloidal fluctuations 
should be rather small for the one-dimensional 
reduced spectra P(/) = P xx + Pyy + Pzz - How¬ 
ever, we note that the influence of poloidal fluc¬ 
tuations increases in the kinetic range and that 
recent findings indicate a slightl y different scal¬ 
ing of the poloidal fluctuat ions ( Chen et al. 2010t 


Podesta fe TenBargel 20121 ). 


As we are able to estimate the scalar function $ 
according to equation 0, we can now proceed to 
describe the diagonal elements of the spectral ten¬ 
sor with a single scalar function T. This function 
depends on the three-dimensional wave vector k 
but since we assume axisymmetry along the mean 
magnetic field, this dependence further reduces to 
'L(k) = fc||). Each wave vector represents a 

fluctuation with a certain wavelength along a cer¬ 
tain direction in space. These fluctuations have 
a constant angle to the mean magnetic field, i.e., 
global and local mean magnetic field are identical 
in our model. For magnetic power spectral densi¬ 
ties, the scalar function T is in units of nT 2 m 3 . 


2.1. Power Spectral Density in Frequency 
Space 

In order to compare modeled energy densities 
in three-dimensional k-space (e.g., equation (U4[) l 
with measured PSD, we need to calculate the 
reduced one-dimensional spectrum in frequency 
space P ( f ) from the spectral tensor given in equa¬ 
tions As the modeled spectra are calcu¬ 

lated with the assumption of statistically homoge¬ 
neous fields and Taylor’s hypothesis, the measured 
data, which the modeled spectra will be compared 
with, need to be statistically stationary. When 
our model is compared to observations, therefore 
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both assumptions, i.e., statistical stationarity and 
Taylor’s hypothesis need to be verified for the par¬ 
ticular properties of the turbulent system and the 
relative velocity of the moving rest frame with re¬ 
spect to the plasma under investigation. For the 
application of Taylor’s hypothesis, it is required 
that k • v w, where k • v is the frequency of 
structures or fluctuations being convected over the 
spacecraft with velocity v and ui is the frequency 
of the corresponding wave or fluctuation in the 
plasma frame. Time-stationarity of observed tur¬ 
bulent fiel ds is often c hecked f or with stationarity 
tests (e.g. lMatthaeus fc Goldsteinl 1982 ). Approx¬ 
imately, it can also be estimated with visual in¬ 
spection of wavelet scalograms. 

In the solar wind, where most of our subsequent 
analysis is applied to, Taylor’s hypothe sis is we ll 
satisfied in the inertial range (Matthacu s fc Goldstein 


fl982i : iNarita et al .11201 .'ll I Howes et al. 2014h . How¬ 
ever, when the Mach number drops below unity, 
e.g., close to the Sun or in planetary magneto¬ 
spheres, Taylor’s hypothesis might not be valid 
anymore, in general. Violation of Taylor’s hy¬ 
pothesis may also happen at kinetic scales, where 
phase speeds of kinetic Alfven waves (KAW) and 
whistler waves can become comparable to the rel¬ 
ative plasma velocity v and frequencies ui in the 
plasma frame strongly increase. 

However, strong anisotropy of the fluctuations, 
e.g., k l fey, may help to satisfy Taylor’s hy¬ 
pothesis: For KAW, the phase velocity is strongly 
reduced for nearly perpendicular propagation as 
WIi.kaw oc 008(0). Therefore, Taylor’s hypoth¬ 
esis holds in the solar wind for a critically bal- 
anced K A W c ascade even in the dissipation range 
( Howes et al.l 12014 ). The same reasoning applies 
to strongly anisotropic KAW fluctuations in the 


magnetospheres of Jupiter an d Saturn (jSaur et al 


2002; von Paoen et al. 2014 1. Our model is also 


applicable for non-relativistic cosmic rays. Here, 
typical energies of lOMeV correspond to ~ 4 • 
lO'ms -1 , for which Taylor’s hypothesis is clearly 
satisfied. Whistler wave turbulence on the other 
hand violates Taylor’s hypothesis and can there¬ 
fore not be described by our model. 

If Taylor’s hypothesis and homogeneity apply, 
the diagonal elements of the reduced spectral ten¬ 
sor in frequency space, Pu(f ), can be obtained by 
integrating the three-dimensional energy densities 
Su( k) over a plane perpendicular to the flow di¬ 


rection v 

oo oo 

Pii(f) =Jdt e l2 * ft j dfc 3 S u (k) e~ ik ' vt 

—oo —oo 

oo 

d 3 k Su(k) 5(2Trf—k x vsin(9)—k z vcos(9)) (5) 

—OO 



(Fredricks fc Coroniti 1976 


Here, we assumed that the z-axis is parallel to 
the background magnetic field (e z || B 0 ) and 
Gy = so that e x lies in the plane spanned 

by B 0 and v, whic h is sometimes call ed the quasi¬ 
parallel direction ( Bieber et al. 1996h . 


In order to evaluate equation (0 numerically, 
we apply a coordinate transformation that reduces 
the dimensions of the integration. If we rotate the 
coordinate system by the angle 90° — 9 about the 
y-axis such that 


K = 

k x sin(0) + k z cos(0) 

(6) 

K = 

k y 

(7) 

k' z = 

—k x cos (9) + k z sin(0) , 

(8) 


then k' x will be aligned with v and the plane of 
integration lies in the y'—z' plane. The unprimed 
parallel and perpendicular wave numbers now take 
the form 


k± = yj[k' x sin(0) — k' z cos(0)] 2 + k'^f (9) 

*11 = k' x cos(0) + k' z sin(0) . (10) 

In the primed coordinate system equation 0 be¬ 
comes 

Pu(f) = l J°° d 3 k' Su{k!)S^-tA . (11) 

Now the Delta function can be evaluated and 
the integration over three dimensions in equa¬ 
tion 0 is thus reduced to two dimensions k' y and 
k' z , which is much easier to calculate numerically. 
Note, that several different wave vectors map to 
the same frequency in spacecraft frame. In fact, 
the power spectral density Pu(f) is determined 
by integration over a plane with normal vector 
k n = {k x sin(0), 0, k z cos(0)). For the numerical 
evaluation, the scalar functions T' and $' are in¬ 
serted in equation m using equations 0-0 in 
the rotated form. With this, we are able to calcu¬ 
late the power spectral densities Pu(f ) in space¬ 
craft frame. 
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3. Results for a Critically Balanced Tur¬ 
bulence Model 


In this section, we introduce a k-space tur¬ 
bulence model ranging from MHD to electron 
scales based on commonly discussed models in 
the literature. With this model we calculate 
one-dimensional reduced PSD for arbitrary mea¬ 
surement geometries. We focus on the theory 
of a strong critically balanc ed tu rbul ent cascade 
as originally proposed by iGoldreich fc Sridliail 
( 19951) for Alfvenic turbulence. It has been fur- 
ther developed to include the ion kinetic r ange 
( Howes et al. 20081 : Schekochihin et al. 2009h an d 


the electron kinetic range (IT enBa rge et al.1120131) . 
Observations indicate that it is suited to de¬ 
scribe solar wind t u rbulence (Horburv et al . 2008; 
Wicks et al. 1201 (it Ichen et al. I 201 ofi ~ although 
these observations have not been explicitly com¬ 
pared with reduced spectra rigorously derived 
from three-dimensional k-space models. However, 
with the forward model derived here, we are able 
to perform such an explicit test and are able to 
check if the observations are truly in agreement 
with critical balance and other suggested proper¬ 
ties of three-dimensional k-space. 


3.1. Critically Balanced Cascade on MHD 
Scales 


According to Goldreicli fe Sridhar ( 19951) . the 
(kinetic) energy distribution on MHD scales can 
be described as 


£(k) 


vl 


k^L 1 / 3 


f U 1/3 - 


, 2/3 


( 12 ) 


where V\ = Bq/^/JIqq is the Alfven speed with g 
the mass density of the plasma, L the outer scale, 
and “/(it) is a positive, symmetric function of 
u, that becomes ne gli gibly small, when |u| /g> 1" 

( Goldreicli fc Sridhar 19951) . It is this function 
/(w), which contains the critical balance relation 
of wave vectors on MHD scales 


L~ 1/3 k 2 / 3 . 


(13) 


The energy is supposed to be isotropically injected 
into the system on the outer scale L , where the ex¬ 
citation is assumed to be strong 5vl ~ V\. This 
relation leads to spectral anisotropy, i.e., differ¬ 
ent scalings depending on the field-to-flow angle 


9. Cho et al. (120021 ) find that the function that 
best explains the magnetic fluctuations in their 
simulation of incompressible MHD turbulence, is 


-Emhd = 0 


Vl 1 /3 


k . 10/3 exp —L 1 / 3 


JM 

U 2/3 


(14) 

According to equation m, this indicates that 
the magnetic fluctuations on the outer scal e L are 
strong, 5B ~ Bq- However, Forman et ah (2011) 


find from fits to solar wind data that Bq overesti¬ 
mates the energy of the turbulent fluctuations by 
a factor of 2 — 3. Therefore, the energy level of the 
modeled power spectra cannot be obtained unam¬ 
biguously. The energy distribution given by equa¬ 
tion m describes the Alfvenic fluctuations and 
we thus set ^ = Lmhd on MHD scales. Accord¬ 
ing to equation GD, this constrains the poloidal 

p 

function to $ = - E r-\I>. 

-r_L 

Next to the exponential function used in equa¬ 
tion ED. there are several different possibilities 
to describe f(u), e.g., Dirac Delta, Heaviside or 
Gauss functions. The particular choice of f(u), 
however, results in similar spectral ani sotropies as 
found in th e case of MHD scale s by Forman et al.l 
( 2011 ) and Turner et al. ( 2012 1) . For 9 = 90°, the 
one-dimensional PSD scales as P(/) oc / _5 / 3 and 
for 9 = 0° as P(f) oc f~ 2 in the inertial range. Ex¬ 
act power-laws in frequency space only occur for 
these two extreme cases. In the literature, how¬ 
ever, it is often implicitly assumed that the spectra 
in the intermediate range 0° < 9 < 90° also follow 
power-laws with spectral slopes 5/3 < k < 2 and 
that the spectral index k of the reduced spectra is 
- similar to the spectral index of three-dimensional 
k-space - independent of frequency. In this paper, 
we will see that these assumptions are incorrect. 


3.2. Transition From MHD to Kinetic 
Scales 

When the observed Doppler-shifted frequencies 
reach scales close to the characteristic ion scales, 
such as the ion gyro radius pi or the ion inertial 
length Ai, the MHD approximation breaks down 
and one needs to take into account kinetic effects 
to describe the plasma dynamics. For the remain¬ 
der of this paper, we use the gyro radius as the 
controlling kinetic scale with an associated criti- 
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cal balance on ion kinetic scales 

fc|, ~ L-^p-^k 1 / 3 (15) 

according to iHowes et al.l (|2008h . Here, we drop 
the factor (ft + 2/(1 + Te/ft)) 1 ^ 6 , where ft is the 
ion plasma beta and T/ and T e the temperatures of 
ions and electrons, respectively, because the factor 
is of order unity in the solar wind and Saturn’s 
magnetosphere. Note, that our forward model can 
similarly be formulated for other theories or with 
the inertial length as the controlling scale. 

We model the transition from MHD to ki¬ 
netic scales as an abrupt change in the function 
£mhd (k) at k±pi = 1 without energy loss, e.g., 
due to ion resonances. Let us denote the func¬ 
tion i?MHD j given in equation Ob as the one 
applicable for kj_pi < 1 and Ekaw as the one 
applicable for k±pi > 1. At k±pi = 1, both 
energy distributions are equal, EMHu(k±Pi=l) = 
£ , KAw(fcj_Pi=l). Also, the kinetic energy distribu- 

_7 /o 

tion, Ekaw, i s assu med to scale as P(f) oc f ± ' 


for 9 = 90° ( Howe s et al. 1200811 . These require¬ 


ments are fulfilled in the expression 


Ekaw = 


B 2 0 


L^p] 


1/3 


u-n/3 

c ± 


exp 


—L 1 / 3 


Pi 


1/3 \k x 


1/3 

a,_l , 
(16) 

this leads to a 


For the parallel cascade {6 = 0 
scaling of P(f) oc f~ 5 . 

At even smaller scales, we approach the elec¬ 
tron gyro radius p e and the energy distribution 
changes once again from Ekaw to an electron or 
dissipation range distribution Fed- Here, Lan¬ 
dau damping is assumed to weaken the cascade so 
that there is no more para ll el transfer of energ y 


1994 Howes et al. 200 


(S ridha r & G oldreich 
The electron or dissipation range fluctuations 
are modeled with an associated critical balance 
201311 of 


(TenBarge et al 


■£- ,/ V v V /a 


(17) 


A functional form that satisfies the equality of 
-Ekaw and -Eed at k±p e = 1 can be given by 


-Eed = 


Bl 


^ 1/3 Pi 


1/3 


k , n ^ 3 exp i — 


(-£ i/s ft / vy 3 ife i) 

(18) 

It can be shown that this leads to a scaling of 
P(f) oc / _8,/3 for the perpendicular cascade (9 = 


90°) and an exponential decay P(f) oc exp(—/) 
for the parallel cascade. 

The energy distribution in k-space according 
to our model is shown in Figure |T] for the pa¬ 
rameters given in Table |Tj It shows logarithmi¬ 
cally equidistant iso-contours of energy densities 
according to equations <m, cud and m in a 
double-logarithmic plot as a function of k x and k z 
normalized by the ion gyro radius pi. The charac¬ 
teristic slopes of the critical balance relations, 2/3, 
1/3, and 0 in the MHD, KAW, and ED regime, 
respectively, can be seen as the boundary, where 
the energy becomes negligible (dark blue). Note, 
that the energy density decreases with k\\ and k ±, 
which leads to a seemingly negative slope in the 
electron dissipation range. 

Recall from equation (0 and subsequent discus¬ 
sion that the reduced one-dimensional spectrum is 
calculated by integrating over a plane given by 

k x sin(0) = — — — k z cos(0) . (19) 

v 

Projections of these planes into the k x —k z plane 
are shown as dashed lines in Figure [l] for an angle 
9 = 1° at logarithmically equidistant frequencies 
/ = 10~ 4 —10Hz. For 9 = 0°, these dashed lines 
would be horizontal and parallel to the ft-axis, 
and for 9 = 90°, they would be vertical and par¬ 
allel to the k z - axis. However, due to the double- 
logarithmic nature of the plot, the dashed lines for 
intermediate angles 9 are curved. Note, that even 
for a very small angle 0 = 1°, the plane of integra¬ 
tion turns quasi-perpendicular to the ft-axis, i.e., 
lies parallel to planes expected from integration 
for 9 = 90°, for sufficiently large k x . 

3.3. Cascade Toward Quasi-Perpendicular 
Spectra 

While equation ([Till can be integrated analyt¬ 
ically for angles 0 = 0° and 9 = 90°, it can only 
be numerically evaluated for intermediate angles 
0° < 9 < 90°. Here, we carry out this numerical 
integration for the first time considering all scales 
from the MHD to the electron kinetic regime. Re¬ 
sults of this calculation using the parameters given 
in Table [T|are shown in Figure [2] We plot the re¬ 
sulting PSD for several field-to-flow angles as a 
function of frequency / in spacecraft frame. Note, 
that frequency can be transformed into normalized 
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-5 -4 -3 -2 -1 0 1 2 3 4 



-5 -4 -3 -2 -1 0 1 2 3 4 


iog(k x p.) 


Fig. 1.— Logarithmically equidistant color-coding 
of energy density in k-space calculated with equa¬ 
tions (fTdl) . (IT6|) . and (fl8| . Thin dashed white 
lines show integration planes after equation (flT)l) 
for logarithmically equidistant frequencies / = 
10~ 4 —10 Hz and a field-to-flow angle 6 = 1°. Due 
to the double-logarithmic plot these planes appear 
as curved lines. Note, that energy resides only in 
wave vectors k > L _1 . 


wave number according to ^-pi = kpt , particu¬ 
larly for parameters in Table [H / a: 0.95 • kpi. 


For 9 = 90° the spectral breaks of the PSD, 
which mark the transition from MHD to KAW 
and KAW to ED range, are located at the Doppler- 
shifted frequencies corresponding to pi and p e (de¬ 
noted by vertical dashed lines). Here, the spectral 
slope of the reduced PSD steepens from 5/3 to 
7/3 and from 7/3 to 8/3, respectively. For smaller 
angles 9 < 90° these spectral breaks are less pro¬ 
nounced and occur at lower frequencies. The ex¬ 
pected spectral slopes for 0 = 0° and 9 = 90° are 
shown as black lines in Figure [2] to guide the eye. 
Here and throughout this paper, we use plasma 
parameters characteristic of the solar wind at 1 AU 
as 


of 

Alexandrova et al. ( 

20091. Schckochihin et al. 

(2009), and Chen et al. 

(201 oh. 


We note a particular interesting feature in Fig¬ 
ure [2j the spectra for oblique angles are steepened 
in a short frequency range around the first spec¬ 
tral break at 2nf/(v sin(9))pi = k±pi ~ 1 and 
then flatten out toward a slope of 7/3 at higher 
frequencies. The smaller the angle 9 , the longer 



Fig. 2.— PSD for different angles 9 as a func¬ 
tion of frequency / with plasma parameters given 
in Table [H The underlying critically balanced k- 
space structure consists of three different regimes 
with spectral indices 5/3, 7/3 and 8/3 for 9 = 
90°. Note, that no damping is included. Verti¬ 
cal dashed lines show Doppler-shifted gyro-radii p; 
and p e for 9 = 90°. At high frequencies all spec¬ 
tra tend toward a slope of 7/3 in KAW and 8/3 
in ED range. Colored dots show the analytically 
approximated transition frequency , where the 
slopes of the spectra approach the slope of the 
perpendicular spectrum. f^ D is calculated with 
equation (IA12D assuming a pure KAW cascade. 


the range, where PSD are steeper than 7/3. Such 
an asymptotic be havior has be e n pre dicted for 
the MHD range bv lForman et al.l ( 201lh . however, 
without specifying the frequency range, where this 
transition could be expected. Here, we show for 
the first time that this transition happens on a 
short frequency range around the spectral break 
in the solar wind. We emphasize the importance 
of this result as it contradicts the commonly as¬ 
sumed constancy of the spectral index with fre¬ 
quency. Also, it implies that PSD in the kinetic 
range are almost exclusively observed with their 
perpendicular slope (here 7/3 and 8/3) if not sub¬ 
ject to damping. In the remainder of this section, 
we explain this feature as a geometrical or sam¬ 
pling effect. 


The flattening of the PSD, although puzzling 
at first, can be explained with the anisotropic dis¬ 
tribution of energy in k-space. Because of the lin¬ 
ear relation between k x and k z in equation (1191) . 
the point of maximum curvature of the integration 
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Table 1: Solar wind parameters used for our model. 


Parameter 

Value 

Background magnetic field B 0 

5nT 

Outer scale L 

10 6 km 

Plasma speed v 

600 kms” 1 

Alfven speed V\ 

60kms _1 

Ion gyro radius pi 

100 km 

Electron gyro radius p e 

2.3 km 


plane in Figure [l] grows faster than the critical bal¬ 
ance relations k\\ ~ in the MHD (a = 2/3), 
KAW (a = 1/3), and ED (a = 0) range. Hence, 
for increasing frequencies /, the plane of integra¬ 
tion given by equation m effectively has only 
contributions from the part quasi-perpendicular 
to the fc^-axis in the double-logarithmic plot of 
Figure [I] Therefore, at frequencies above a tran¬ 
sition frequency fi d, the power spectrum at a 
given field-to-flow angle 8 will have a slope within 
a small error Ak to that at 0 = 90° appropriate 
to MHD (5/3), KAW (7/3), or ED (8/3). The 
transition frequency //d thus marks the bound¬ 
ary, where the anisotropy of the k-space energy 
distribution becomes so large that it is indistin¬ 
guishable from a 2D-distribution in the reduced 
PSD within measurement error Aft. It is this fea¬ 
ture that causes the spectral slope of the PSD to 
asymptotically approach its perpendicular value 
at high frequencies. 

A derivation of an analytically approximated 
expression /| D (see equation (IA12D 1 can be found 
in Appendix [A] The frequencies, where the PSD 
turn quasi-perpendicular, are shown in Figure[2]as 
colored dots assuming a pure KAW cascade. Note, 
that the transition frequency // D , in contrast to 
the spectral break, shifts to higher frequencies for 
smaller angles 8. Assuming a plain MHD cas¬ 
cade, the approximate transition frequency from 
equation (1A12I) gives /2 d,mhd(^ = 20°) = 0.2 Hz. 
The corresponding transition frequencies for 8 = 
[1°, 5°] reach into the kinetic range of scales. The 
transition to quasi-perpendicular scaling is always 
found at frequencies below the electron spectral 
break, k±p e ~ 1, which is why frequencies are 
not shown in the ED range. 

We find that the transition for KAW toward a 
quasi-perpendicular cascade for almost all field- 


to-flow angles (namely 8 >4°, see Figure Qj] 
in the appendix) occurs already below 1 Hz, 
where the spectral break is observed at typi¬ 
cal solar wind conditio ns ( Sahraoui et al. 2010; 
Alexandrova et all 1 2012 1. We conclude that spec¬ 


tra of a critically balanced KAW cascade without 
damping on kinetic scales will almost exclusively 
be observed with spectral indices close to 7/3 and 
8/3 in the ion and electron kinetic range, respec¬ 
tively. Only for small angles and in a short fre¬ 
quency range between the first spectral break and 
1 Hz a measurably steeper slope might be observed 
in the solar wind. 

We stress the importance of this result: we 
have shown that the spectral index in a criti¬ 
cally balanced cascade for oblique field-to-flow 
angles 0° < 8 < 90° is not constant with growing 
frequency. Instead, it evolves toward a quasi¬ 
perpendicular slope as a function of frequency. 
Therefore, we cannot expect to observe signifi¬ 
cantly steeper slopes than 7/3 and 8/3 over a 
broad frequency range in the KAW and ED range, 
respectively, for a critically balanced cascade in 
the solar wind. Indeed, we see in Section 13.41 
that such slopes only appear if the fluctuations 
are additionally subject to damping. This transi¬ 
tion of reduced PSD toward quasi-perpendicular 
spectra contrasts a common understanding in the 
literature, where it is often implicitly assumed 
that the PSD have a power-law shape with a 
constant spectral index that only depends on 
the field-to-flow angle. It is interesting to note 
that a three-dimensional direct num eric al simula- 
tion with strong guide field by Grappin & Muller 
( 2010 ) produced a qualitatively similar transi¬ 
tion toward a perpendicular cascade with a 8- 
independent slope, although the flattening of their 
spectra does not follow equation (IA12I) . 


3.4. Anisotropic Damping 


Several mechanisms have been proposed to 
explain the dissipation of fluctuations on small 
scales, e.g., ion-cyclotron dampin g, Landau damp- 


ine and current sheet formation (Matthaeus et al. 

1990; 

Leamon et al.l 1999; Dmitruk et al. 2004; 

Howes 

2009; Schekochihin et al. 2009; TenBaree et al. 

2013). 

Most of the proposed mechanisms are 


anisotropic with regards to the background mag¬ 
netic field. In this section, we analyze how possible 
anisotropic damping terms that act on the energy 
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Fig. 3.— Ratio of PSD without damping (P 0 ) 
and with ion-cyclotron damping (Pj c ) according 
to equation (BUI) as a function of frequency / for 
several field-to-flow angles 9. Colored dots show 
/f D according to a pure KAW cascade. 


distribution in three-dimensional k-space change 
the characteristics of reduced one-dimensional 
spectra. Here, we assume that damping leads 
to an exponential decay of the energy in k-space 
associated with the damped wave vectors. 


From our assumption of critical balance, it 
follows that k± fc||. Therefore, damping at 
a fixed scale affects the turbulent cascade first 
through perpendicular wave vectors and only 
later through wave vectors parallel to the mean 
magnetic field. It has been argued that the 
ion-cyclotron resonance at fc|| ~ D lc /Vx 1 where 
f2ic = — is the ion-cyclotron frequency, has only 
minor influence on a critically balanced cascade 
as it is reached only at very high perpendicular 
wave numb ers k \ , where Landau damping already 
domin ates (iHowes et al. 2008| Schekochihin et al 


20091 : Cranmer fc van Ballegooiien 2012 ). With 


our synthetic spectra, we are now able to test this 
argument quantitatively. 


3-4-1- Damping of Parallel Wave Vectors 


For the ion-cyclotron resonance, the damping 
rate is low for k\\V a -C D- lc and i ncreases rapid! 
as k \ \V\/D.\ c . approaches unity ( Howes et al.l 200 
Cranmer fc van Ballegooiienl l2012f) . We may 


therefore model the ion-cyclotron resonance as 


an exponential decay according to 

P ic (k)=P(k)-exp(-^^) . (20) 

To quantify the effective damping, we model 
damped spectra Pi c using the parameters in Table 
U| and find only minor changes compared to their 
undamped counterparts Co¬ 
in Figure |H we show the ratio Po/Pi c as a func¬ 
tion of normalized wave number for several an¬ 
gles 9. Both Pq and Pi c cover all three ranges 
from MHD to KAW to electron kinetic scales. 
For quasi-parallel spectra (9 = 1°), a peak fac¬ 
tor of 1.34 is observed around the spectral break. 
Power spectra Pi c for larger field-to-flow angles 
9 > 5° only differ a factor of 1.2 from PSD with¬ 
out damping. On electron scales, the difference 
between damped and undamped spectra remains 
constant as there is no more parallel transfer along 
the cascade. The asymptotic factor Po/-Pic for 
2irf/V > pf 1 and 9 = 90° can be calculated by 
integration of equation (fl8l) over k z according to 

^2nf/V>pf 1 ,9=90°) 

MC 

/ d k z exp (-P 1 / 3 p l 1/3 pe /3 |fc z |') 

— OO 

OO , x. 

/ d k z exp (-A 1 / 3 p l 1/3 p e /3 \k z \- \k z \^) 

—oo ' ' 

« 1.21 . ( 21 ) 


This shows that the influence of cyclotron damp¬ 
ing on the spectral shape can be neglected in the 
solar wind for critically balanced turbulence. 


However, one might think of systems, in which 
the ion-cyclotron resonance does change the form 
of the PSD significantly. Such systems could be 
characterized by (1) low plasma density, which 
leads to large Alfven speeds; (2) a critical balance 
exponent a close to unity; and/or (3) a much 
smaller outer scale than the L ~ 10 9 —10 lo m 
found in the solar _ wind ( Howes et al.l 2008; 


Schekochihin et al. 20091) . so that high fc|| values 
are reached earlier in the cascade. An enhanced 
ion-cyclotron damping then leads to a visible de¬ 
crease of spectral power around the spectral break. 


3-4-2. Damping of Perpendicular Wave Vectors 


Let us now t urn to da mping of perpendicu¬ 
lar wave vectors. Alexandrova et al. (|2012) found 
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that reduced power spectra measured in the solar 
wind can empirically be described by 


P(k±) oc k ± K exp(— k±p e ) (22) 


with a spectral index of n = 8/3. This result is 
consistent with numerical gyrokinetic sim ulations 
( Howes et al. 2011 : TenBaree et all [201 .'ll) and we 
may treat the exponential damping term of equa¬ 
tion (1221) as a proxy for electron Landau damping. 
Note, that we apply the damping term of equa¬ 
tion m to three-dimensional k-space according 
to 

-Edamp(k) = E( k) • exp(—fcj_p e ) . (23) 


The corresponding distribution in k-space is 
shown in Figure [U where both damping of par¬ 
allel and perpendicular wave vectors (equations 
and m) is included. While the power at a 
certain frequency is obtained by integration over a 
two-dimensional plane, the damping term depends 
on k± and is therefore axisymmetric with respect 
to k z . This means that the integration along k y 
includes perpendicular wave numbers larger than 
k± = 27r//(usin(0)). Consequently, the effect of 
the exponential damping term is stronger when 
applied to three-dimensional k-space compared to 
a reduced one-dimensional spectrum. Therefore, 
we are able to produce similar results to those of 
equation (l22l) with k = 8/3, although we apply 
a less steep spectral index of only 7/3 in the ion 
kinetic range (see Figure [5]). 


3-4-3. Reduced PSD Subject to Damping 

Figure [5] shows damped PSD according to equa¬ 
tions and (1231) with the same plasma parame¬ 
ters (TableQ]) as in Figured As expected, the high 
frequency part is strongly affected and the damp¬ 
ing leads to a characteristic exponential decay. 
The spectra are dominated by damping of per¬ 
pendicular wave vectors and the damping terms 
are found to be generally more effective at small 
angles 6. The latter can be understood consider¬ 
ing the plane of integration given by equation m- 
The power at a certain frequency primarily stems 
from perpendicular wave numbers k±. For small 
angles 6 , this involves larger k±, which are more 
strongly damped than a corresponding spectrum 
with angles close to 90°. 

Although strong damping sets in not before 
electron scales in wave vector space, the KAW 



Fig. 4.— Exponentially damped and critically bal¬ 
anced energy distribution in k-space according to 
equations and (l23l) . The distribution con¬ 
tains three ranges for MHD, ion and electron ki¬ 
netic scales. Thin dashed white lines show inte¬ 
gration planes after equation 03 for logarithmi¬ 
cally equidistant frequencies / = 10~ 4 —10 Hz and 
a field-to-flow angle 0 = 1°. 


range of the PSD, p” 1 < k < pj 1 , is already af¬ 
fected: Figure [5] shows the undamped spectrum 
for 9 = 90° as black dashed line and it is fairly 
visible that the spectrum subject to damping is 
steeper. In fact, we measure a spectral slope of 
2.63 (shown by black dashed line) in the ion ki¬ 
netic range, 1 < kpi < 42 for Pdamp(0=9O°), al¬ 
though the corresponding energy in k-space scales 
with 7/3. 

If T.i/T e decreases or the electron gyro radius 
increases, the exponential decay moves to lower 
frequencies and further steepens the spectra. The 
two opposed mechanisms, flattening toward per¬ 
pendicular slope and damping, can cancel each 
other and may lead to a range of seemingly con¬ 
stant slope. This shows that an approximate 
power-law of the measured PSD is not equivalent 
to the absence of damping. Further, it shows that 
the observed spectral slope of the reduced spec¬ 
trum on ion kinetic scales is not necessarily the 
spectral index predicted for this range by the un¬ 
derlying theory. 
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Fig. 5.- PSD subject to damping according to 
equations KU1) and (0^)1 as a function of frequency 
/ for parameters given in Table [Q The black 
dashed line shows the PSD for 9 = 90° without 
damping. Colored dots show for orientation at 
what frequency the spectrum would turn quasi¬ 
perpendicular for a pure KAW cascade without 
damping. 


4. Application to Solar Wind Observa¬ 
tions 

By numerically evaluating equation (ED, it 
is possible to calculate power spectral densities 
Pu(f) for any component i = x,y,z and given 
plasma parameters (plasma speed, outer scale, 
held-to-flow angle, background magnetic field, 
Alfven speed, gyro radii). Here, we compare 
our results to in-situ measur e ment s in the solar 
wind mad e by Horbur v et a.1.1 ( 20081) in the MHD 
range an d Chen et al. (2010 ) in the kin etic range. 
Podesta ( 20091) and Wicks et al. ( 2010h have pre¬ 
sented similar results in the MHD range, w h ich ar e 


consistent with those from iHorburv et al.l ([2008). 


We show that the measured spectral anisotropies 
are in accordance with our model and can thus be 
described by a critically balanced cascade. How¬ 
ever, the influence of damping turns out to be 
more important than previously thought. 


4.1. MHD Turbulence 

From the assumption of critical balance it fol¬ 
lows that the reduced PSD scales as P a f ~ 5 / 3 
for 9 = 90° and P oc f ~ 2 for 0 = 0°. Fig¬ 
ure [G] shows exemplary for the case of MHD, i.e., 
equation ED only, that the transition between 



Fig. 6.— Spectral anisotropy in the inertial range 
for different values of the outer scale L for MHD 
scaling, i.e., integration of ED over expression 
ED- The spectral index is obtained in frequency 
range / = [10~ 4 —0.1] Hz. A smaller outer scale 
results in a transition of the spectral index from 2 
to 5/3 at larger angles 9. 


these two scalings is controlled by the outer scale 
L. Here, we have calculated the spectral slope 
from a least-squares fit in the frequency range 
/ = [10 -4 , 0.1] Hz. Note, that the spectral index 
slightly varies in this frequency range and, there¬ 
fore, no exact power-law is observed. However, the 
result shown in Figure [6] impressibly characterizes 
the influence of the outer scale on the cascade. 
It can be seen that the perpendicular scaling is 
reached only at very large angles 9 for small outer 
scales L. This reflects the evolution of the turbu¬ 
lent cascade along the critical balance path. While 
the energy is isotropic at the outer scale, it grows 
increasingly anisotropic as it cascades to smaller 
scales. 


We can use this functional dependence of the 
spectral anisotropy to check whether our critically 
balanced k-spac e distribution g e nerat es results in 
agreement with Horburv et al. ( 20081) and, if so, 
which value of L fits the observations best. For 
several outer scales L, we calculate the spectral 
anisotropy and fit the spectral indices k(L,9) to 
their results. To evaluate the critically balanced 
power spectra, we use the complete set of energy 
distributions, -Emhd, Ekaw and -Eed, together 
with anisotropic damping, i.e., exponential damp¬ 
ing at the ion-cyclotron resonance and at electron 
scales according to equations (f20l) and (l23l) . re¬ 
spectively. 

[Horburv et al.l ( 20081) used magnetic field data 
from the Ulysses spacecraft at 1.4 AU in the fast 
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solar w ind (V — 750 km s - 1 ). From lMcComas et al 
(l2000l) . we estimate the proton gyro radius as 
Pi = 190 km and due to lack of a better estimate, 
we assume T) = T e to calculate the electron gyro 
radius. For comparison, we also calculate the best 
fit for a slab+2-D turbulence model, where we use 
a spectral index of 2 for th e slab component an d 
5/3 for the 2-D component (jHorburv et al.l 1201 111 . 
To quantify the goodness of fit we use a reduced 
error 


X = 


\ 


1 N 

-Y 

N ^ 

i=1 


(Hi ~ Ki ) 2 


(24) 


where k,; and eh are the spectral indices and their 
corresp onding errors , respectively, taken from 
Horburv ct al. ( 2008lf . N is the number of angle 
bins, and Ki the modeled spectral indices, which 
are obtained by a least squares fit in the same 
frequency rang e, 15—98 mHz, as those used by 
Horburv et al. ( 20081) . The spectral anisotropies 
in the frequency range 15—98 mHz for the best fit 
parameters are shown in Figure [ 7 ] 

In general, the fit is better for the critically 
balanced turbulence model (% = 2.8—3.0) than 
for the slab+2-D turbulence model (x = 4.6). 
This can be seen in Figure [7] where we show the 
observed and modeled spectral anisotropies. We 
find that an outer scale of L = 10 9 m in the un¬ 
damped case and L = 10 10 m for the damped cas¬ 
cade give the best results. These values are in ac¬ 
cordance with observations of the transition from 
/ _1 -spectra to / _5 / 3 -spectra in the solar wind, 
which is belie ved to mark the end of the energy in¬ 
jection scale ( Schekochihin et al. 2009). The spec¬ 
tra at small angles 9 are found to be steeper than 
k = 2. This is not caused by damping. The differ¬ 
ence between damped and undamped cascade is, 
according to Figures [Hand 01 very small. Instead, 
the steep slopes result from the fact that the fit¬ 
ting range / = 15—98 mHz includes the steeper ki¬ 
netic range cascade, where slopes are significantly 
steeper than k = 2. For 9 = 5°, e.g., / = 98mHz 
corresponds to k±pi = 1.8. 

For slab+2-D turbulence, we find the best fit 
has 30% slab and 70% 2-D turbulence, which is 
also in ac cordance with resu lts obtained in the so¬ 
lar wind ( Bieber et al. 19961) . Note, however, that 
we use two different spectral indices for slab and 
2-D turbulence while Bieber et al. 0199 6) used the 
same slope for both slab and 2-D turbulence. 



Fig. _ 7.— Spect ral anisotropy found by 

Horburv et al. ( 2008i ) (black crosses), as obtained 


from our model with an outer scale of L = 10 9 m 
for undamped (black line) and L = 10 10 m for 
damped (red line) critical balance, as well as for 
slab+2-D turbulence with 30% slab (blue line). 


The anisotropy of the power, P(9)/P(9= 5°), 
measured at a fixed frequency of / = 61 mHz is 
shown in Figure 01 The results of the damped 
critically balanced cascade well reprod uce the 
anisotropy found bv lHorburv et al. ( 20081) for an¬ 
gles 9 < 40° but overestimate the power at larger 
angles. The undamped cascade generally shows 
a weaker anisotropy but is in agreement with 
the measurements at 9 ~ 90°. Similar to what 
has been found for the spectral index, a critically 
balanced cascade fits the data much better than 
slab+2-D turbulence. 

Given the simplicity of our p arameters, the 


presen ted fits to the results of iHorburv et al 


( 2008 ) can be regarded as a qualitatively success- 
ful reproducti o n of t he observed data. Although 
Forman et al.l ( 201ll ) showed that the observed 
spectral anisotropy is in agreement with a critical 
balance on MHD scales, this is the first time that 
these results have been analyzed with a model in¬ 
cluding kinetic range scalings and damping terms. 
The quality of the fits i ndicates that th e spect ral 
anisotropy observed by [Horburv et al.l (120081) is 
consistent with a critically balanced cascade. 


4.2. Kinetic Range Turbulence 

For an an alysis of kine tic ran ge magnetic field 
fluctuations, I Chen et ah ( 2010 ) used data from 
CLUSTER during fast solar wind conditions with 
moderate ion plasma /3j ~ 1. They calculated the 
spectral index using a wavelet based method de- 
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Field-to-flow angle 0 


Fig. 8.— Power spectral anisotropy 

P(8)/P(8= 5°) at frequency / = 61 mHz. 

The results of damped and und amped critical 


balance are close to the results oflHorburv et al 


(120081 1 while a combined model of slab and 2-D 
turbulence underestimates the anisotropy. 


scribed by iHorburv et al. ( 2008ll . Ion and elec¬ 
tron gyro radii as well as th e bu lk plasma ve - 
locity are given in Table 1 of Chen et al . 1 d2010h 
and the power anisotropy is given as P±/P\\ ~ 20. 
The spectral anisotropy is calculated separately 
for parallel and perpendicular fluctuations. Here, 
we limit our fit to the perpendicular fluctuations 
P±. The spectral anisotropy of parallel fluctu¬ 
ations can not be reproduced su fficient ly well 


which is consistent with findings of (Forman et al 
(120131 1. For the modeled spectra, we choose for 
the outer scales L = 10 9 m for the undamped and 
L = 10 10 m for the damped ca scade based on our 
fit to the IHorburv et al' ( 2008 1 data. 

In Figure [9l we show the results obtained from 
our for ward calc ulation and those presented in 
Idien et al.l ( 2010h for the ion kinetic range 1.5 < 
kpi < 6. Additional to the damped KAW cas¬ 
cade with ft = 7/3, we show results for a damped 
critically balanced cascade with an assumed spec¬ 
tral index of 8/3 and critical balance exponent 
a = 1/3. This scenario is included as hypotheti¬ 
cal case for visual orientation only. Although the 
modeled results for damped and undamped spec¬ 
tra do not fit the data within error tolerances, 
the damped cascades show qualitatively similar 
spectral anisotropies. Here, the variation from 
0 = 0° to 8 = 90° is smooth and not as abrupt 
as in the undamped case. For cross-comparison, 



Field-to-flow angle 0 


Fig. 9. — Spectral anisotropy of fluctuations per¬ 
pendicular to the local magnetic field derived by 
Chen et HI (l201di ) (black crosses) for kinetic scales 
kpi = [1.5,6] and numerically determined results 
for critically balanced KAW cascades. The spec¬ 
tral anisotropies for the damped cascades (red line 
for ft = 7/3 and blue line for k = 8/3) fit better 
to the measured data compared to the undamped 
cascade (black line). 


we point the reader to Figure [5] where it can 
be seen that PSD are steeper for smaller field-to- 
flow angles 8. The undamped cascade, in con¬ 
trast, leads to spectral slopes of ~7/3 at almost 
all angles (cf. Figure [2]). Even at 8 = 5°, where 
the observed spectral index is ~3.3, the spectral 
index of the undamped cascade is already 2.4. 
This is the consequence of the transition toward a 
quasi-perpendicular slope as was elaborated in sec¬ 
tion 13.31 The spectral index of the damped KAW 
cascade on the other hand is 3.2 at 9 = 5°, and 
thus much closer to the observation. 


The results in Figure [9] show that an undamped 
cascade leads to a sharp increase of the spectral 
index at small field-to-flow angles, which appears 
rather step-like. This result has interesting impli¬ 
cations for the interpretation of critically balanced 
plasma turbulence. It is often assumed that a crit¬ 
ically balanced cascade in the ion kinetic range 
can lead to any spectral slope, 7/3 < ft < 5, for 
intermediate field-to-flow angles 9. However, we 
have shown in Section 13.31 and in equation (IA12I) 
that the spectral index is ^7/3 for most angles. A 
smooth or slow variation of the spectr al index, as 
seen in the results of Chen et al . ( 20ld ), can there¬ 
fore not be caused by critical balance alone. In- 
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stead, the effect of damping is essential to obtain 
steeper spectra at non-zero angles, which means 
that the observed spectral anisotropy is to a large 
degree determined by the damping mechanisms. 


5. Modeling Spectral Densities at Saturn 


Recently, it has been shown that magnetic fluc¬ 
tuations in S aturn’s plasma sheet f orm a turbu¬ 


lent cascade ( von Papen et al. 2014 b These fluc¬ 


tuations have different properties compared to the 
solar wind. They are embedded within Saturn’s 
strong background magnetic field and propagate 
as Alfven waves along magnetic field lines un¬ 
til they are reflected by density gradients at the 
plasma sheet or the ionosphere. Measurements by 
the Cassini spacecraft during its first seven orbits 
around Saturn indicate that observed power-law 
spectra can be interpreted as a critically balanced 
cascade of KAW. On MHD scales, however, power 
spectra indicating an Alfven wave cascade with 
Kolmogorov-like power-law are observed only spo¬ 
radically because large-scale magnetospheric pro¬ 
cesses dominate this range of scales. As Taylor’s 
hypothesis holds in Saturn’s magnetosphere for ki¬ 
netic Alfven waves with k± >> , we are able to 

test the observations with our forward model and 
to investigate the spatial distribution of observed 
spectral indices as a function of radial distance to 
Saturn. 

As the turbulent cascade is predominantly ob¬ 
served on kinetic scales, we restrict our model to 
the reproduction of kinetic range spectra. We an¬ 
alyze the radial distri bution of spectral indice s re¬ 
ported in Figure 12 of lvon Panen et ah ( 2014f ) and 
focus on an explanation for the change of spectral 
slopes inside 9 R S: where 1 R s = 60268 km is the 
planetary radius of Saturn. Inside 9 i? s , increas¬ 
ingly flatter spectra are observed, which so far 
could not be explained. As there are two distinct 
electron populations in Saturn’s magnetosphere, 
we also analyze which dissipation scales control 
the onset of damping. Comparison of modeled and 
observed spectra may thus shed more light on the 
physics of the turbulence and may help to under¬ 
stand which electron population is energized by 
the turbulent cascade in Saturn’s magnetosphere. 

To model synthetic spectra in Saturn’s magne¬ 
tosphere, we u se the parameters listed in Table [2] 
as reported in Ivon Papen et al .i ( 2014 b where we 


use 2 Hw ~ L as outer scale. These parameters 
fluctuate in time and vary strongly with radial dis¬ 
tance to Saturn. The basic plasma parameters at 
Saturn - namely velocity, scale height, ion tem¬ 
per ature and densi t y - ar e based on observations 
by iThomsen et al.l ( 201Clh for water group ions, 
which are the main constituent of Saturn’s mag¬ 
netosphere. To calculate the electron gyro radii, 
which control the onset of the empirical damping 
term in equation (123[) , we use the electr o n tem pcr- 
ature models given bv ISchippers et al~ ( 20081 Ta- 
ble 1) derived from CAPS/MIMI measu rements 
( Young et, al.l l2004 iKrimigis et al. 2004 b The 
field-to-flow angle 9 is calculated assuming the 
plasma flow to be in azimuthal direction only. We 
further use the measured power anisotropy (-&-), 


averaged in the kinetic range, to estimate the ratio 
T/d* between the energy of toroidal and poloidal 
fluctuations according to equation |4|). The slopes 
of the modeled spectra are calculated from the 
trace of the synthetic spectral tensor. 


Table 2: List of model input parameters to calcu¬ 
late synthetic spectra in Saturn’s magnetosphere. 


Parameter 

Description 

B 

Magnetic field 

9 

Field-to-flow angle 

V 

Plasma velocity 

H w 

Scale height of water group ions 

pw 

Water group ion gyro radius 

Pe 

Electron gyro radius 

V A 

Alfven velocity 

(ft) 

Power anisotropy 


Equatorial tempera tures have been derived by 


ISchippers et al. ( 2008h for the cold or thermal 


(<100 eV) and hot (>100 eV) electron popula¬ 
tions. These temperatures control the electron 
gyro radii according to 


Pe,c/h 


y/2m e k B T etC / h 

eB 


(25) 


where T ec and T e h denote the cold and hot electron 
temperatures, respectively. The radial profiles of 
the temperatures and densities are shown in the 
top and bottom panels of Figure flOl respectively, 
together with the water group ion temperatures 
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Fig. 10. — Equatorial electron temperatures (top) 
and densities (bottom) in Saturn’s magnetosphere 
for the hot (red) and cold (blue) electron popu¬ 
lations after [Schippers^etj^ (2003) ■ For compari¬ 
son, we also show the equatorial temperatures and 
densi ties of water group io ns (black) as obtained 
from Thomsen et al.1 ( 2010h . 


and densities for comparison. The model for the 
hot electrons is valid for distances up to 18i? s . 
The temperature of the cold population is more 
variable outside 15i? s and thus only provided in¬ 
side of that distance. However, due to the lack of a 
better estimate, we approximate the cold electron 
populations also outside of 15 R s with this model. 
Both temperatures similarly peak at around 9 R s . 

In the following, we model reduced spectra us¬ 
ing electron gyro radii from the cold and hot elec¬ 
tron populations and compare both to the ob¬ 
served spectra. The larger electron gyro radius 
of the hot electron population leads to an earlier 
onset of damping, i.e., the damping term is larger 
and its effects are already important at lower fre¬ 
quencies compared to damping on cold electron 
scales. This causes the slope of the hot electron 
spectrum to be slightly steeper than the cold elec¬ 
tron spectrum. 


In order to compare our modeled results with 
the observations, we need to add noise to the 
synthetic data and assume that the excitation of 
turbulence is not strong, i.e., we assume Svl < 
Va- The latter is as expected from turbulence 
excited within a stron g planetary magnetic field 
( von Papen et al EoTi . The noise levels are com¬ 
puted as the sum of magnetometer instrument 
noise rq, quantization noise n qi and aliasing. The 
instrument noise is measured and must therefore 
be added to the spectrum before calculating the 
aliasing. This leads to a synthetic spectrum P' 
given by 

P\f) = P(f)/cE + rn/f + n q 

+A(f)-[P(f)/c E +n i /f} , (26) 


where P(f) is the synthetic spectrum without 
noise, is a factor, which considers the weaker 
excitation discussed above, and the aliasing func¬ 
tion 


m = £ 


(///»y)“ 


71=1 


[2n - (///ny)] 


7 ^ 

n — 1 


(///ny)* 


(///ny)] K 

(27) 

is determined according to IPodesta et al. (2006f) 
for a power-law of k = 7/3 and Nyquist frequency 
/ ny . In our case, this function takes on a value 
of 1.3 at / = / ny = 3.6 Hz and drops rapidly for 
lower frequencies. 

The second term on the right-hand-side of equa¬ 
tion m is the instrum ent noise, w hich is nj = 
75pT 2 /Hz accordi ng to Dougherty et al. ( 20041) 
and goes with 1// (Russell 19721 ). The third term 
is the quantization noise, which can be estimated 
as n q = ^AB 2 At ( Russelll[l972 1 with sampling pe¬ 
riod At = 0.14 s and the quantization of the mag¬ 
netic field data, e.g., AB — 4.9 pT in FGM range 0 
(IDoughertv et a.l.ll2004ll . The factor ce effectively 
corrects the energies of the synthetic spectra to 
the observed ones, however, it does not change 
the spectral index. The factor ce = ( Pobs/P ) is 
on average 50 and is calculated as the geometri¬ 
cal mean of the ratio of observed P 0 bs to synthetic 
spectra P in the normalized wave number range 
k±pw = [2,50], where k ± = and p w is 

the water group ion gyro radius, and for signal-to- 
noise ratios SNR > 5. 

We carry out a forward modeling for each of 
the 1136 observed 10 min time series and compute 
the spectral slopes in the same frequency ranges 
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as determined for the observed data, i.e., in the 
range 2 < k±pw < 50, and for SNR > 5 of th e 


respective measurement ( von Paoen et al. 2014f ). 


As an example, we show an observed spectrum 
from 9.5 R s and its corresponding synthetic spec¬ 
tra in Figure [IT] The black line shows the ob¬ 
served spectrum, the blue line the spectrum with 
dissipation controlled by the gyro radius of the 
cold electrons and the red line the spectrum with 
dissipation controlled by the gyro radius of the hot 
electrons. All of these spectra include noise with 
exception of the red dashed line, which shows the 
spectrum P° ot without noise for comparison. 

Inside the fitting range, depicted by two vertical 
dashed lines, the three spectra are very similar. 
Their spectral indices are given in the legend of 
the figure and the modeled results (fv co id = 2-32 
and Khot = 2.61) are close to the observed slope 
of k = 2.46. However, at high frequencies the 
spectrum controlled by the gyro radius of the hot 
electrons clearly fits the observation better. At 
frequencies higher than the fitting range, the noise 
level shown by a dotted line leads to a flattening 
of the synthetic spectra. On average, this leads 
to a slight decrease of the spectral index within 
the fitting range by 0.06 and 0.07 for cold and hot 
electron spectra, respectively, compared to spectra 
without noise (see red dashed line in Figure 1X11) . 

Figure [12] shows the radial distribution of ob¬ 
served spectral indices in Saturn ’s magnetosphere 
as black crosses (see Figure 12 in von Papen et al 


(2014j)). The spectral slopes for damping on cold 
electron scales (blue dots) are generally less steep 
than the observations and close to the undamped 
spectral index of 7/3. This shows that the corre¬ 
sponding synthetic spectra have not yet reached 
dissipation scales, where the spectral energy de¬ 
creases significantly. Damping on hot electron 
scales, on the other hand, leads to steeper spec¬ 
tra that agree much better with the observations. 
Clearly, the decrease of the electron temperature 
inside 9 R s leads to shallower spectra because the 
damping is reduced for smaller electron gyro radii. 
This indicates that the change of slopes inside 9 R s 
can be explained by damping effects. However, 
the radial profiles of both electron temperatures 
are nearly identical so that the difference between 
the two gyro radii can in principle be compen¬ 
sated by a simple factor c > 1 in the exponent 
of the damping term, i.e., exp(—c- k±p ec ) in equa¬ 



Fig. 11. Power spectral density of a 10 min time 
series measured in Saturn’s magnetosphere at a 
radial distance of 9.5 R s (black line). The blue 
and red lines show synthetic spectra for damping 
controlled by the gyro radius of cold and hot elec¬ 
trons, respectively. Associated PSD are labeled 
-Fcoid and Phot, respectively. Vertical dashed lines 
show the fitting range in which the spectral index 
is calculated. The determined values are given in 
the legend. The dotted line shows the noise level 
corresponding to Phot and the dashed red line the 
spectrum P® ot without noise. 


tion (EH) . Note, that for larger distances (> 12 R s ) 
the spectral indices Khot become shallower than 
the observations, which indicates either the pres¬ 
ence of additional effects on the turbulent cascade 
or magnetospheric variations. 

In summary, the forward modeling of kinetic 
range spectra has provided a possible explanation 
for the change of spectral slopes inside 9 R Sl which 
is solely based on first principles of a KAW cas¬ 
cade and an empirical term to describe damping 
on electron scales. Our results indicate that the 
energy transferred along the kinetic range cascade 
is preferably deposited into the hot electron pop¬ 
ulation. The heating examined here could thus 
be the primary process, which maintains the hot 
electron population in Saturn’s magnetosphere. 

6. Conclusion 

We present results from a numerical forward 
model to evaluate reduced PSD from from an 
arbitrarily distributed energy density in three- 
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Distance [Fy 


Fig. 12.— Radial distribution of observed spectral 
indices in Saturn’s magnetosphere shown as black 
crosses. Spectral indices of synthetic spectra are 
shown as dots for damping on cold (blue) and hot 
(red) electron scales. Damping on hot electron 
scales explains qualitatively the change of slopes 
inside of 9R S . 


dimensional k-spacc0. Given a critically balanced 
k-space distribution of energy, we investigate the 
functional dependence of the reduced spectra on 
several parameters, such as the field-to-flow an¬ 
gle or the outer scale. Such an analysis has been 
carried out for the first time covering the com¬ 
plete range from MHD to electron kinetic scales. 
We show that for intermediate field-to-flow an¬ 
gles 0 , the spectral slope of undamped critically 
balanced turbulence is not constant, i.e., the re¬ 
duced PSD is not an exact power-law (see Fig¬ 
ure E2). Instead, the reduced spectra evolve with 
frequency in the undamped case toward a quasi¬ 
perpendicular spectrum. This is a pure sampling 
effect of highly Doppler-shifted measurements at 
any angle 9 > 0° caused by the anisotropy k± 
k || increasing with frequency. Power spectra in this 
quasi-perpendicular cascade range have a spectral 
slope corresponding to the spectral index of the 
perpendicular cascade fc(d=90°). For PSD that 
are additionally subject to damping (Figure 0, 
this transition to quasi-perpendicular spectra is 
masked by Landau damping on electron kinetic 
scales. Here, the spectra become steeper and 
steeper at higher frequencies. 

The transition frequency, where the change of 
slope to k k k( 0=9O°) is reached for the un- 


2 The forward code is available in the sup- 
plementary online material and under 
www.geomet.uni-koeln.de/en/research/turbulence 


damped cascade, can be approximated by a sim- 
pie expression f% D given in equation (IA12I) . Under 
typical solar wind conditions, /f^ for the ion ki¬ 
netic regime is smaller than or on the same order of 
the observed spectral break at /b ~ 0.3 Hz. There¬ 
fore, significantly steeper slopes than the perpen¬ 
dicular spectra on kinetic scales can only be ex¬ 
plained by additional damping effects. We find 
that an empirical dampin g ter m at e le ctro n scales 
of the form exp(— k±p e ) (Alexandrova et al. 120121 ) 
already affects the spectral index in the ion kinetic 
range p< k < p~ l of the reduced spectrum and 
steepens it measurably. 

We apply our model to in-situ measurements 
in the solar wind and show that tu rbulent fluctu¬ 
ations measured by iHorburv et al. 


(1200811 


are m 


good agreement with a critically balanced cascade 
and less so with a combined slab+2-D turbulence. 
In the kineti c range we f ind good agreement with 
the results of lChen et al.l ( 2010lf , which can be ex¬ 
plained by a damped critically balanced KAW cas¬ 
cade. While the spectral break and damping terms 
have only m inor influence on the fit to the ob¬ 
servations of IHorburv et al. ( 20081) . the inclusion 
of a damping term is essential to explain the re¬ 
sults in the kinetic range. Indeed, spectral slopes 
for an unda mped cascade differ strongly from the 
IChen et all ( 2010 1 results and are approximately 
7/3 for nearly all field-to-flow angles in the ion 
kinetic regime. This is caused by the transition 
toward a quasi-perpendicular slope. This means 
that damping is an integral part of the kinetic 
range cascade and the dominating factor in the 
observed spectral anisotropy. 

The forward modeling technique is further ap¬ 
plied to observations of turbulent magnetic field 
fluctuations in Saturn’s magnetosphere. Here, we 
show that the measured reduced spectra are qual¬ 
itatively in agreement with a critically balanced 
kinetic Alfven wave cascade, w hich furt her cor¬ 
robora tes the interpretation of Ivon Papen et al~ 
2014T) . The observation of shallower spectra in¬ 
side 9 R s can be reproduced by damping on elec¬ 
tron scales controlled by the hot electron popula¬ 
tion. This indicates that the dissipation of tur¬ 
bulent magnetic fluctuations predominantly heats 
the hot electron population at Saturn. This addi¬ 
tional insight on magnetic turbulence in Saturn’s 
magnetosphere shows the advantage of applying a 
first principles forward modeling technique to de- 
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A. Estimation of Transition Frequency to Quasi-Perpendicular Spectra 


Here, we derive an approximate expression for the frequency /f^, at which the reduced power spectrum 
of an undamped critically balanced cascade turns quasi-perpendicular, i.e., the PSD is characterized by a 
spectral slope k = n(9= 90°) + An with a small error Ak. For mathematical tractability we derive this 
approximate expression under additional assumptions. Later, we show by comparison with our full three- 
dimensional model that the resultant expression still provides a good estimator when these assumptions 
are relaxed. The assumptions we use for our derivation are a two-dimensional k-space (fc||, k±) and a 
critically balanced cascade that is controlled by a Dirac Delta function instead of the exponential function 
i n equ ation ©■ The latter seems to be a strong simplification but iForman et alJ (1201 lh and Turner et a l 


( 2012h have shown that there are only minor differences between the spectral anisotropies of PSD controlled 
by an exponential and a Dirac Delta function. 

In order to estimate the implications of our assumption of a two-dimensional k-space, it is important to 
clarify the three-dimensional geometry of the problem. In Figured] we see a slice along the k x -k z plane. 
While the planes of integration (white dashed lines) are two-dimensional and expand straight in k y direction, 
the energy distribution is three-dimensional and axisymmetric with respect to k z . This means that the plane 
of integration intersects with both Dirac Delta surfaces (for k z > 0 and k z < 0) in a complicated curve, 
which depends on outer scale L , field-to-flow angle 9 , and critical balance exponent a. However, most 
of the energy along these curves resides at k y = 0 and decreases with growing wave number in y-direction. 
Therefore, we assume that the assumption of a two-dimensional wave vector space is suitable for our following 
derivation. This is subsequently corroborated by the agreement of the such derived expression with our fully 
three-dimensional model results. 

In two dimensions, the energy distribution of critical balance controlled by a Dirac Delta function may 
be written as 

E 2D (k ± ,k ll )(xkl K 6(\k ll \-L- 1 / 3 p a - 2 / 3 k^j , (Al) 

where k = [5/3, 7/3, 8/3] are the one-dimensional spectral indice s and a = [2/3 ,1/3,0] the critical balance 
exponents on MHD, KAW, and ED scales, respectively (see also Grappin fc Miillei 2010t h For the sake of 
simplicity, we use only one kinetic scale p in equation (EH): in the KAW range p = p z , while the result for 
the ED range can be obtained by setting p = pip e . The distribution of energy density in k-space given by 
E 2 D is shown in Figure [13] The intersection of the two Dirac Delta branches with the line of integration 
given by equation dD are two points: 


ifcii = L~ 1/3 p a - 2/3 k' 


^_L > 


(A2) 


positive fcy i(/cx i i) and negative fcy 2 (^, 2 )- The locations of these points are depicted as red circles in Figure 
□31 Field-to-flow angles smaller than 90° yield £q_,i < fc_L,o < fcj_, 2 , where 


k_ L,o 


2t Tf 
v sin(0) 


(A3) 


is the wave number corresponding to k\\ = 0. Note, that we do not need to make any assumptions regarding 
the energy at zero parallel wave number. The two locations fcj_,i /2 are uniquely defined by equations (fl9l) 
and (1A2I) . It follows that 


fc_L,l/2 


2tt/ 
vsin(0) 
2 7T/ 
vsin(0) 


-F ^||,i /2 cot(0) 

=F L- 1/3 p a - 2/3 /c“ 1/2 cot(0) . 


(A4) 


With the particular choices of a = [2/3,1/3], this equation can be written as a cubic polynomial with real 
solutions fc_L,i/2(0)- Even though a cubic algebraic equation can be solved analytically, we search for an 
approximate solution to achieve a mathematically simpler expression, which is more easily to work with. 
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E(k^,k N ), 6=1 




'°9io< k x> 


Fig. 13.- Double-logarithmic plot of energy distribution given by a Delta correlated critical balance model 
as given in equation ED (thick black line) in k-space with blue line showing the integration curve for a 
field-to-flow angle 0 = 1°. Red circles indicate intersections of integration curve with Dirac Delta branch and 
thin black lines denote the locations of k_ l,i/ 2 - The blue circles show the locations approximated according 
to equation (1X71) . Dashed lines show the locations of k±o and kug. 


Let us now look at the energy E 2 d in k-space that contributes to the integration: 


E 2 £> oc fej_^ + i 


(AS) 


which is essentially the summation over both critical balance branches. In the following, we show that this 
energy asymptotically approaches Eq oc 2with increasing frequency, thus leading to a scaling according 
to 9 = 90°. At a certain frequency /f^, the measurement uncertainty will be larger than the difference 
between the energies E 2 d and Eq and, therefore, the observed PSD will scale like the perpendicular cascade. 
We proceed to estimate the frequency for which the reduced PSD turns quasi-perpendicular by demanding 
that the ratio E 2 d/Eq be almost unity. Because the spectral index is obtained from the logarithms of the 
power spectral energy, we demand 


In 


-i- k~ K 
ft X,l + *'±,2 

Z/C _L,0 


< e ) 


(A6) 


where e > 0 is chosen according to typical measurement uncertainties Aft for spectral slopes in solar wind 
observations. 

To proceed further, we approximate 




fc ||0 = T- 1 /3 p -V3 fc 


a 

_L,0 • 


(A7) 


The geometrical interpretation of this approximation is shown in Figure 1131 Instead of the exact solutions 
shown by the red circles, we now evaluate the integral at the approximated locations shown by the blue 
circles. These are slightly shifted along the integration curve toward lower fc_L-values. Now, we can write 
equation (|A4I) as 

k±, 1/2 « 1l,o T T _1/3 p a_2/3 fe“ i0 cot(0) 

= fc_L,o =F A . (A8) 
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Inserting this in equation (IA6I) . we get 


1 - 


kj l,c 


+ 1 


k_ L,C 


< 2e £ 


(A9) 


If we expand the two terms on the left-hand-side of equation (IA9I) to second order around A/k±p = 0, we 
find that 

A 


k±, t 


< 


2e e — 2 
k(k ± 1 ) 


(A10) 


Inserting A = L 1 / 3 p“ 2//3 fc“ 0 cot(0) and k±p = 27 if / (usin(0)) into this equation finally yields 


. usin(6M i 2-3a ( 2e £ — 2 

/ > —-- 


2n 


c(k + 1 ) 


cot(0)i- 0 


(All) 


For an appropriately chosen e, the difference between energies Eid and Eq for all frequencies / > /| D is 
so small that the scaling will be quasi-perpendicular within measurement errors An. Note, that the change 
of slope in frequency space is a pure sampling effect and does not mean that the nature of the turbulent 
cascade in wave vector space is changing at this frequency. For e <C 1, we can further simplify equation (IA11I) 
to determine the approximate transition frequency that marks the boundary to quasi-perpendicular scaling 
as 


/Id 


v sin(0) 
27T 



2e 


k(k ± 1 ) 




(A12) 


Note, that equation (IA12I) can additionally be used to test whether an observed steep slope in a reduced spec¬ 
trum of critically balanced turbulence can be explained by spectral anisotropy only, i.e., without damping, 
or if additional damping is needed to explain the steep slope. 

In order to check in which turbulent range lies, i.e., in the MHD, the KAW or the ED range, the 
following procedure might be used: One chooses the appropriate critical balance exponent for the observed 
range of scales, namely a = 2/3 for MHD, a = 1/3 for KAW, or a = 0 for electron kinetic scales. Inserting 
a in equation (IA12I) together with the relative plasma velocity v, observation angle 9, outer scale L, and 
controlling kinetic scale p (p = pi for KAW and p = y/pip e for ED), one obtains the transition frequency 
/fb fo r the observation geometry under consideration. All frequencies / > /lb will be characterized by a 
quasi-perpendicular spectrum. 

To estimate in which range (MHD, KAW, ED) the transition toward quasi-perpendicular spectra occurs, 
one can calculate the corresponding perpendicular wave number according to k± = 2nf2 D /(v sin(0)). If the 
resultant k± falls outside the range used to calculate f% D (with certain a and p), then the transition toward 
quasi-perpendicular spectra will either happen on smaller scales (if k± is too large for the corresponding 
range) or the transition has already happened on larger scales (if k± is too small). The latter case, e.g., is 
found for the ED range in the solar wind as the transition toward quasi-perpendicular spectra already occurs 
on KAW scales. To find the correct transition frequency, the procedure then needs to be repeated with the 
a and p values from the alternative range. 

Note, that the error parameter e must be determined according to the measurement accuracy, which 
usually requires modeling the corresponding spectra. We find that e = 0.05 in equation (IA12I) represents an 
uncertainty of the spectral index of An = ±0.1 in the KAW and Ak = ±0.02 in the MHD range for typical 
solar wind conditions. For the KAW range, this uncertainty is on the order of the measurement error in 
solar wind observations. 

In the following, we test the approximate expression given in equation (1A12I) against our three- 
dimensional model results. For that matter, we determine the transition frequencies from modeled PSD 
without a break at k±pi ~ 1 or k±p e ~ 1, i.e., for equations m and (USD separately, each extending over 
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Fig. 14.— Transition frequencies as a function of 9 for critical balance values a = [2/3,1/3, 0]. Solid lines 
show calculated from the approximate expression in equation (IA12I) using e = 0.05 and crosses show 
f 2 D d obtained numerically from modeled PSD for an uncertainty in the spectral index of Are = 0.02 on 
MHD and Ak = 0.1 on ion kinetic scales. 


the whole range of wave vectors. Because of divergence at low wave number kj_, this method is numerically 
not stable if we use an energy distribution in k-space, which is solely given by equation (1181) . From the 
modeled PSD in the MHD and KAW range, we calculate the spectral slopes between two consecutive data 
points (fn+i/fn ~ 1-3) asK = re(0=90°) + Are, where re(0=90°) = [5/3, 7/3] in the MHD and KAW regime, 
respectively. We then define f£$ d as the frequency, where Are decreases below thresholds of Are = ±0.1 in 
the KAW and Are = ±0.02 in the MHD range. 

In Figure [TH we show the such derived f^ d (crosses) for parameters given in Table |T] compared to the 
analytical approximation f% D (solid lines) given in equation (IA12I) for critical balance values a = [2/3,1/3, 0] 
as a function of field-to-flow angle 0° < 6 < 90°. The modeled three-dimensional results agree well with the 
estimations from equation (1A12I) for e = 0.05, which are shown by solid lines for the MHD, KAW, and ED 
range. On MHD scales, the transition toward quasi-perpendicular spectra occurs over a broad frequency 
range and is therefore not as apparent as on kinetic scales. 

The calculated transition frequencies f% D are shown in the corresponding figures throughout this paper 
as colored dots for the KAW range. The transition frequencies f% D MHD for the MHD regime according to 
equation (IA12I) are not shown in the figures as they usually fall outside the MHD range, i.e., f% D MHD > ^r- 
However, for 9 = 20° equation (IA12D results in f? lD MHD = 0.2 Hz. The transition frequencies according to 
the ED range always appear at lower frequencies than the electron spectral break and are therefore not 
shown. 
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